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1 Introduction 


The existence and properties of two-dimensional conformal field theories crucially depend 
on the value of the central charge c, which is a parameter of the Virasoro symmetry 
algebra. For example minimal models, the only consistent theories whose spectrums are 
made of finitely many irreducible representations, exist for the values 


c = 1 — 6 


(p - qf 

pq 


( 1 . 1 ) 


where p > q > 2 are integers. Liouville theory, the simplest nontrivial theory with a 
continuous spectrum, is universally believed to exist for c > 25. Moreover, it is widely 
believed that Liouville theory has a consistent analytic continuation to all values of c 
except c < 1. (See [1] for a review.) 
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Liouville theory (continuation) 


25 


Liouville theory 


Minimal models 


For c < 1, the analytic continuation of Liouville theory does not have a well-dehned 
limit. On the other hand, the crossing symmetry equations for the three-point function 
behave similarly for c < 1 and for c > 25: the degenerate equations, that is those of the 
crossing symmetry equations which can be written explicitly, have a unique solution in 
these cases. For c > 25, that unique solution is the Liouville three-point function [2]. For 
c < 1, the unique solution is known [3-5], but it is not known whether it solves the rest 
of the crossing symmetry equations, and therefore corresponds to a consistent theory. 

Recent hndings on the two-dimensional Potts model provide renewed motivation for 
studying Liouville theory with c < 1. In its formulation in terms of random clusters, 
the Q-states Potts model comes with a parameter 1 < Q < 4. Special cases include 
percolation {Q = 1) and the Ising model {Q = 2). In the critical limit, the clusters 
become conformally invariant random fractals, and the parameter Q is related to the 
central charge c with 0 < c < 1. (See [6] for a review.) While the critical exponents 
of the Q-states Potts model have been known for a long time, not much is known on 
the correlation functions. Recently, it was observed that probabilities for three points to 
belong to the same cluster are given by certain values of the three-point function which 
solves the degenerate crossing symmetry equations [7-9]. This suggests that this three- 
point function should belong to a consistent theory. Our aim is to propose a precise 
dehnition of that theory. 

The crossing symmetry equations involve not only the three-point structure function, 
but also the spectrum of the theory. So we should determine the spectrum of Liouville 
theory with c < 1. It was suggested that there exists a timelike Liouville theory [10], where 
conformal dimensions are bounded from above. We will dismiss this as inconsistent, and 
focus on the hypothesis that conformal dimensions are bounded from below, as is the 
case for the other values of c. With this hypothesis, correlation functions appear to be 
singular, but the singularities can be eliminated by a simple regularization. 

The crucial evidence that our ideas are correct will be provided by numerical checks of 
crossing symmetry. For that evidence to be acceptable, the corresponding computer code 
should be publicly available. One could even argue that the software which is used should 
be free and open source. So we wrote Python code in the IPython Notebook format, and 
published it on GitHub. In particular, the code for generating our data tables and plots 
can be found in the notebook “article_support. ipynb”, which has the same structure 
as the present article. 


2 



Using our code, we checked crossing symmetry not only in Liouville theory with c < 1, 
but also in Liouville theory with generic complex values of c, and in generalized minimal 
models. We also checked crossing symmetry in certain non-rational theories which exist 


at rational values of c (1-1), and which generalize the c = 1 Runkel-Watts theory. And 


we tested a few speculations of our own on the existence of further theories at rational 
values of c. 

2 Liouville theory in the conformal bootstrap approach 

Let us sketch the conformal bootstrap approach to Liouville theory, and explain how the 
results depend on the value of c. More details and references can be found in the review 
article [1]. 

2.1 Definition and spectrum of Liouville theory 

We dehne Liouville theory as a two-dimensional conformal held theory which obeys two 
main assumptions: 

1. The spectrum contains continuously many Verma modules of the Virasoro algebra, 
each of them having multiplicity one. 

2. The theory depends smoothly on the central charge c, moreover correlation functions 
depend smoothly on the conformal dimensions of the helds. 

Let us spell out the meaning and hrst consequences of these assumptions. The spectrum 
is a representation of the symmetry algebra QJ x which is made of two copies of the 
Virasoro algebra T1 with the same central charge c - these copies are called left-moving 
or holomorphic for QJ, and right-moving or anti-holomorphic for ill. The spectrum is 
built from Verma modules Va of the Virasoro algebra, which depend on the conformal 
dimension A. Our assumptions imply that the spectrum is diagonal. 



( 2 . 1 ) 


We still have to specify the allowed values of the conformal dimension A. For c > 1, these 


values can be guessed from some features of the representation theory of the Virasoro 
algebra, which we now review. 

The set of representations of the Virasoro algebra at a given value of c is a ring for 


an associative product called the fusion product, which is useful in conformal held theory 


because it constrains the operator product expansions of the corresponding helds. Verma 
modules have the rather trivial fusion product 



( 2 . 2 ) 


However, there exist degenerate representations {'JZ{r,s)}r,seN*, whose fusion products in¬ 
volve only hnitely many terms. These fusion products have simple expressions provided 
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we parametrize the central charge c in terms of the background charge Q and coupling 
constant b, 


c — 1 + 6Q^ , Q — 6 + 6 ^ , (2-3) 

and the conformal dimension A in terms of the momentum a, 

A = a{Q-a) . (2.4) 

The fusion product of the degenerate representation 7l(r,s) with a Verma module is a sum 
of rs Verma modules, 


r —1 ,8 — 1 

2 2 


T^{r,s) X Vq = ^2 ^2 




■ 1 — r • 1 — s 

2 2 


where the sums run in steps of 1. Algebraically, a degenerate representation is the quotient 
of a reducible Verma module by a nontrivial submodule, 

K,.,., = . (2.6) 

^( —r,s) 

where we introduce 

. (2.7) 

To conclude on Virasoro representations, let us mention that unitary Verma modules exist 
only if c > 1. More precisely, if c > 1 then Va is unitary if and only if A > 0. 

For c > 1, we will use the assumption of unitarity for guiding us in guessing the 
spectrum. A hrst guess would be to include all Verma modules which obey the unitarity 
bound A > 0. However, we have also learned that the natural variable of Virasoro 
representation theory is not the conformal dimension A, but the momentum a, and our 
assumption of smoothness of correlation functions really means analyticity in a. In terms 
of a, the unitarity bound becomes 


A>0 a e (0,(5) U ( Y+ 


It would be difficult for correlation functions to be analytic if a was allowed to live in the 
segment (0, Q) with its hnite endpoints. So we exclude this segment from the spectrum, 
and the allowed values of the momentum become 


c- 1 


a e Y + * . 


We will later see that this is true not only for c > 1, but for all c G C. For the moment 
let us introduce a real variable P for the momentum, which is related to A and a by 


Q 

a = j + iP 


( 2 . 10 ) 
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2.2 Crossing symmetry equations for correlation functions 

Solving Liouville theory means computing its correlation functions. An A^-point correla¬ 
tion function is associated to N fields, where a field is an operator-valued function on the 
Euclidean plane. There is a state-field correspondence between states in the spectrum, 
and fields. In a Verma module Va, all states are obtained by acting with creation oper¬ 
ators on the module’s primary state. We call Va (or equivalently Va) the primary field 
which corresponds to that primary state. Then the Wpoint function of primary fields 


N 

i=l 



( 2 . 11 ) 


is a function of the central charge, the momentums and the complex coordinates Zi 
on the Euclidean plane. The fields which are obtained by acting with creation operators 
are called descendent fields, and their correlation functions are deduced from correlation 
functions of primary fields through Ward identities. 

In principle, correlation functions only involve fields which correspond to states in the 
spectrum, whose momentums take values as in eq. (2.9). However, it is fruitful to consider 
fields which correspond to more general states, and in particular to states in degenerate 
representations. So we assume that there exists a degenerate primary field V)^,*) which 
corresponds to the primary state of the degenerate representation 7l(r,s)- (We are not 
assuming that V(r,s) can be obtained from Va by analytic continuation, see the discussion 
in Section 3.1.) The existence of a vanishing descendent in the representation 7Z(r,s) (2.6) 
implies the existence of a vanishing descendent field for V(r,s)- This in turn implies that a 
correlation function involving V(^r,s} obeys a Belavin-Polyakov-Zamolodchikov differential 
equation of order rs, whose solutions correspond to the rs terms in the fusion product 
T^{r,s) X Va (2.5). 

The main axiom which underlies the calculation of correlation functions is the exis¬ 
tence of an associative operator product expansion. 


Vai (^l) Vq.2 ('^ 2 ) 



da C{ai, 0 ( 2 , Q 




\cVa{z2) . 


( 2 . 12 ) 


This expresses the product of two fields Va^{zi)Va 2 {z 2 ) as a linear combination of the 
primary field Va (obtained for £ = £ = 1 ) and its descendent fields ££14, where 
£ is a creation operator in the universal enveloping algebra of the left-moving Vira- 
soro algebra. The coefficients in this linear combination are factorized into the univer¬ 
sal factors S'fj q, 2 |q,(^ 1 ) ^ 2 ); which are determined by conformal symmetry and normal¬ 
ized so that Qa^a 10 (^ 15 ^ 2 ) = 1 , and the three-point structure constant C'(ai,Q! 2 ,Q — a). 


^cxi,(y2\<^ 

Here and in the following, the notation 

modulus squared, but for a product of left-moving anc 
2 


^ai,a2|a(^li^2) 


does not stand for the complex 
right-moving quantities, namely 


0)1 ,Oi2 | q ; 


(^ 1 , 2 : 2 ) 


= ^fl,a2|a(^li^2)^4,„2|a(^l42). 

Repeatedly using the operator product expansion, any correlation function can be 
written as a combination of three-point structure constants, and universal quantities. For 
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example, a four-point correlation function can be written as 



/ dag C[ai, a2, Q — as)C{as, a^, a/j^) 





(2.13) 


where we introduced the four-point, s-channel conformal blocks ~ combina¬ 
tions of and similar universal quantities, summed over C. Such conformal 

blocks are determined by conformal symmetry and therefore in principle known, just like 
characters of representations. The calculation of conformal blocks is discussed in Section 
A.2. Here we only need the following properties: 

• J^a\ai\zi) is a meromorphic function of b and a, whose 6 -poles occur for the “min¬ 
imal” values 6 ^ G Q<o (IT), and whose a-poles occur for “degenerate” values 

ex CX(r,s) (2.7). 


• For large values of A = a(Q — a), conformal blocks behave as 

J^i^\ai\zi) ~ (16g)^ with |g| < 1 , (2.14) 

|A|^oo 

where g is a function of 2 : 1 , • • • ^ 4 . 


Our decomposition of a four-point correlation function was obtained using the opera¬ 
tor product expansion of but we could alternatively use the expansion of lAjlAg, 

and obtain the t-channel decomposition. 


( 




dat 0 ( 0 : 2 , 03 , Q — at)C{at, CX 45 cti) 



(2.15) 


The equality of the s- and t-channel decompositions is called crossing symmetry, and 
is equivalent to the associativity of the operator product expansion. (In particular, no 
further constraints arise from the u-channel decomposition [11].) Crossing symmetry can 
be represented using the following graphs, where each vertex corresponds to a three-point 
structure constant: 


/ 


dus 




(2.16) 


Crossing symmetry provides an unwieldy set of equations for the three-point structure 
constant C(q:i, 0 : 2 , exs): we have as many equations as choices of values of cxi, • • •CX4, each 
equation is quadratic, and has inhnitely many terms. Simpler equations, with hnitely 
many terms, can be obtained if we replace one held, say with a degenerate held 
V(r,s). Then cx^ and at take rs possible values, which are determined by the fusion rules 
(2.5). These values are in general outside y -f-iM, which is not a problem since correlation 
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functions are analytic in a. The simplest nontrivial examples are obtained for (r, s) G 
{(1, 2), (2,1)}. With V( 2 ,i), we obtain the following equation: 


E 





(2.17) 


This leads to relations involving C(Q:i,Q: 2 ,tt 3 ± |) and C{ai ± ^,a 2 ,a 3 ), which actually 
imply 

C{ai,a2,a3 + b) _ b~‘^'''^'y{2ba3)'y{2ba3 + 

C(q;i, 0(2, tts) 7(60(3 ± 6 (q;i — y) ± 6(a2 — §)) 

where the function 7 (x) is dehned from Euler’s Gamma function by 


7 ( 1 ) 


r(x) 

r(i-x) ' 


(2.19) 


If we had used the degenerate field V(i^ 2 ) instead of V( 2 ,i), the equation we would have 
obtained would be related to eq. (2.18) by 6 —)■ b~^. Summing up, we have two equations 
of the type 


G(q(i, 0(2, as + b) 
C{ai, 0 ( 2 , 0 / 3 ) 


known , 


G(o(i, 0(2, as + 6 
G(ai, 0 ( 2 , as) 


known , 


( 2 . 20 ) 


which we call degenerate crossing symmetry equations. (These equations appear to 
be linear, whereas crossing symmetry equations should be quadratic in the three-point 
structure constant. Actually, degenerate values of the structure constant of the type of 
G((l, 2), 0 (, 0 ( ± |) are hidden in the right-hand sides.) 

2.3 Solving the crossing symmetry equations 

The degenerate crossing symmetry equations determine how the three-point structure 
constant behaves under shifts of momentums by b and b~^. Let us represent b and b~^ as 
vectors in the complex plane: 


i 


0 1 


b eiR 
c< 1 



b G C 
CG C 




( 2 . 21 ) 


6 G M 
c > 25 


Here and in the following, 6 G C should be understood as 6 ^ M U zR, and similarly 
c G C should be understood as c ^ (—cxo, 1) U (25, cx)). For a generic value of 6 G R U iR, 
the vectors b and b~^ are aligned and incommensurable. Therefore, there exists a unique 
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smooth solution of the shift equations, up to an overall momentum-independent factor. 
We will call this solution after Dorn-Otto and Zamolodchikov-Zamolodchikov if 

c > 25, and if c < 1. For b ^ C, the vectors b and b~^ are not aligned, and there 
exist non-constant smooth functions which are periodic in both b and b~^. Multiplying a 
solution with such smooth functions yields more solutions. 

A three-point structure constant must solve not only the degenerate crossing symme¬ 
try equations, but also the non-degenerate equations. If c < 1 or c > 25, uniqueness of 
the solution of the degenerate equations suggests that it has good chances of solving all 
equations. Moreover, for c > 25, Liouville theory is believed to exist and to be consistent, 
based in particular on semi-classical approaches near the c —)■ -|-oo limit. On the other 
hand, for c G C, it would be very implausible that each one of the inhnitely many solutions 
of the degenerate equations solves all crossing symmetry equations. The existence of more 
than one solution of these vastly overdetermined equations would already be implausible. 

We will now describe the three-point constant for c > 25, and explain how it 

also provides a solution of Liouville theory for c G C. To solve the degenerate crossing 
symmetry equations, we manifestly need functions whose shifts by both b and b~^ generate 
Gamma functions as in eq. (2.18). And indeed there exists a function Tb(a;) such that 


Tb{x + b) ..^bx /u ^ A U--1 /X'. 

—^1 ^^^r^{bx) and ^ = b ^ ^ 7 (f) 


( 2 . 22 ) 


Tb(a:) ^ Tfe(a:) 

This function can be dehned either by an integral formula, which is valid for 0 < Kx < ^Q, 


\ogTb{x) 



sinh^ — x) t 

sinh bt sinh i 
0 


(2.23) 


or by a product formula, which is valid for all x G C, 

Q 


= \ 


n / 

m,n=0 


X 


with f{x) = (1 - x^)e^ 


Y -h mb + nb~^ 


(2.24) 


where Aj, is an unimportant constant. The product formula makes it clear that Tb(x) is 
analytic on C, with inhnitely many simple zeros: 


Tf,(x) = 0 77 X G {-bN - b-^N) u{Q + bN + . (2.25) 


And the integral formula allows us to compute the large momentum behaviour 

logTfed + zP) = -PHog\P\ + lp^ + o{P^) . (2.26) 

^ P—^oo / 

The function T;, is well-dehned except for b G iM, where the product and integral formulas 
diverge. 

Using the function Tb(x), we now write the three-point structure constant in terms 
of the variable P (2.10), 

nliU(Q+2iP.) 

n±,± u (I+*(u ±P2± fs))' 

- 8 - 


^DOZZ(p^,p^^p^) 


(2.27) 









The symmetry of P 2 , P 3 ) under permutations of {Pi} follows from the identity 

T,(x) = T,(Q-x). 

Like the function this three-point structure constant is well-dehned for all c except 
c < 1. A four-point function can then be computed using its s-channel decomposition 
(2.13). Let us check that the integral over G M converges. To begin with, the poles of 
the integrand are safely away from the integration line. There are indeed poles from the 
conformal blocks at 

Pse\J±l{Q + bN + b-^N) , (2.28) 

± ^ 

and poles from the structure constants at 

Pse [j (^±Pi±P 2 ±*(§ + 6 M + riM)) 

±,±,± ^ 

IJ IJ ( ± P 3 ± P 4 ± z (f + 6 N + b-^N) ) . (2.29) 
±,±,± ^ 

Drawing a discrete set of poles P -|- i{bN + 6 “^N) as a wedge, 



the poles of the integrand form 9 wedges on each side of the integration line, leaving a 
strip of width free of poles: 


—Pi — P 2 + P 4 P 3 — Pi —Pi — P 2 0 Pi — P 2 —P 3 + Pi —P 3 — Pi Pi + P 2 

- 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - > 

Ps 


It remains to be checked that the four-point function, as an integral over P*, converges at 
Pg = ± 00 . Using the behaviour of the conformal blocks (2.14), and that of the three-point 
structure constants, which follows from eq. (2.26), we hnd that the integrand behaves as 

2 




Pa 


Pp\Pi\Zi 


rsj 

Pfi—^dlCO 


kl 


2Pf 


(2.30) 
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with |g| < 1 . So the integral converges. This means that the four-point function is well- 
dehned, as computed from its s-channel decomposition. The same holds for the t-channel 
decomposition. Assuming that these two decompositions agree for c < 25, they must 
agree for all values of c except c < 1 , because of the 6 -independent convergent behaviour 
of the integrals, and the fact that the integrands are analytic in b. In other words, the 
DOZZ three-point structure constant obeys crossing symmetry wherever it exists. (See 
also the discussion in [ 12 ].) 

We checked this numerically for various values of the parameters. For example, with 
{Pi) = (1.3,1.01, 0.45, 0.22) and {zi) = (0.26, 0, oo, 1), we hnd the following values for the 
real parts of the s- and t-channel decompositions, depending on the central charge c: 


c 

s-channel 

f-channel 

precision 

36.74 

256250.53844 

256245.59396 

1.9 X 10-5 

17.55 

1104.88358 

1104.85666 

2.4 X 10-5 

2.12 

2.04664 

2.04659 

2.6 X 10-5 

3.00 + 4.00Z 

-2.97769 

-2.97757 

3.6 X 10-5 

0.50 + 2.00i 

-5.92002 

-5.91968 

4.3 X 10-5 

-1.00 +2.lOi 

-11.43323 

-11.43175 

8.3 X 10-5 


Details such as the normalizations of the four-point functions can be found in the code. 
What matters is the agreement between both channels, up to a small relative difference 
which we call the precision. That difference can be made smaller by increasing the nu¬ 
merical accuracy of the computations, see the Appendix. 


3 The problem with c < 1 

Having reviewed the dehnition and solution of Liouville theory for generic complex values 
of the central charge, we now turn to the case c < 1 . 


3.1 The three-point structure constant 

In the case c < 1, the degenerate crossing symmetry equations (2.20) have the unique 
smooth solution [3-5] 


C^^\PuP2,P3) = 




nti ri,a + 2F.) 


(3.1) 


where we introduce the parameter 


/3 = G M . (3.2) 

This solution has sometimes been rejected as unphysical [3], because in the limit Pi —)■ 
iy = |(/d~^ — /3) the resulting three-point function does not vanish for P 2 7 ^ ±^ 3 , and 
therefore cannot be interpreted as a two-point function. This means that the held limVk 
differs from the degenerate held V(i^i) oc identity, which corresponds to the primary state 
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of the irreducible coset (2.6). Instead, lirnV^ corresponds to the primary state 

of the reducible Verma module So we have a non-degenerate held with zero 

conformal dimension, which however has nothing to do with logarithmic conformal held 
theory, since the Virasoro generator Lq is diagonalizable in the reducible Verma module 
Vo (as in all Verma modules). This contrasts with Liouville theory with generic c, where 
limVo, happens to be degenerate, due to the analytic properties of the DOZZ three-point 
structure constant. Nevertheless, as already pointed out in [12], there is no reason to 
expect that limVo, should always be degenerate. 

The structure constant is a meromorphic function of the momentums Pj, with 
countably many poles. In particular, the purely imaginary values of the momentums do 
not play any distinguished role, and there is nothing intrinsically timelike about 
although it has sometimes been called the timelike Liouville three-point function [9]. 
Our notation instead emphasizes the domain where it is the unique solution of the 
degenerate crossing symmetry equations. 

The structure constant obeys the following equation, which involves two different 
values of the parameter (3, 


16- 


2^C<1 

' ^(v^/3) 


\/2P, 


v/ 2 ’ 


^ 2 ) 


rtc<i 

^ 0 ) 


P — 


p \ p<c<l I _p J_ _ J. 

« / (/ 3 ) \ 4 / 3 ! df 


4 / 3 ’ 4/3 




{V2(P + f + j;)) 

(^/2(F + i)) 


(3.3) 


where u/s is a function of (3. This relation is similar to a relation for the DOZZ structure 
constant, which is an important ingredient in proving modular invariance of Liouville 
theory on the torus [13]. The other ingredient of the proof is a relation for conformal 
blocks, which is true for all values of the central charge by analyticity of the conformal 
blocks. The proof then relates modular invariance on the torus to crossing symmetry on 
the sphere - the structure constants in eq. (3.3) are relevant for a torus one-point function 
(numerator) and sphere four-point function (denominator): 


\/2P- 

(\/ 2 /d) 



(3.4) 


So, if we prove crossing symmetry of Liouville theory with c < 1, it will follow that the 
theory is consistent not only on the sphere, but also on the torus and therefore on all 
oriented Riemann surfaces. 


3.2 Why timelike Liouville theory does not exist 

Let us discuss the spectrum of Liouville theory with c < 1. In contrast to the case c > 1, 
there is no continuum of unitary representations from which to build the spectrum. So 
we do not have any guidance from unitarity. 
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However, we already know the three-point structure constant and conformal 

blocks. In the s-channel decomposition of the four-point function, 


( 



[ dPsC'^^\P^,P2,Ps)C<^-\-Ps,P3,P4) 

J E 




(3.5) 


the only unknown ingredient is the domain E where the momentum Pg should be inte¬ 
grated. Let us discuss what this domain can be, based on the properties of the integrand. 
The integrand is a locally analytic function of Pg. It has poles from the conformal blocks, 
which are given by the general formula (2.28). For c < 1, that formula amounts to 

Pse\J±^ (/3(1 + N) - r'(1 + N)) . (3.6) 


If 7 ^ Q, then these poles are dense in the real line. There are also poles from the 
structure constants C^-^, which in contrast to the poles (2.29) from do not depend 

on the values of Pi, • • ■ P 4 . Instead, these poles are again on the real line, 

Pse[j±^[f3N + (3-^N) -{0}. (3.7) 


So all the poles of the integrand are on the real line. But the real line is where we 
should integrate Pg, if we assumed that the spectrum of Liouville theory was given for 
c < 1 by the same expression (2.9) as for the rest of the values of c. It was proposed that 
for c < 1 there exists a timelike Liouville theory instead [10], where the momentum Pg 
would take imaginary values: 




0 


This proposal however has to be discarded, because the integral in the four-point function 
(3.5) diverges at Pg = ±ioo. Using the behaviour of the conformal blocks (2.14), and that 
of the three-point structure constants, which follows from eq. (2.26), we indeed hud 


C"=^i(Pi, P2, P,)C"=^'(-P„ P3, P 


s 


{Pi\Zi] 


\2Pf 


Ps—fitioo 


(3.8) 


which is formally the same behaviour as for generic values of c with Pg —)■ ±cx) (2.30). So, 
in timelike Liouville theory, it is impossible to compute four-point functions. 

One might worry that our argument condemns all timelike theories. After all, the 
divergent factor is universal, as it comes from the conformal blocks. However, the 
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integral expression for the four-point function diverges only if the s-channel momentum 
Pg is allowed to explore large imaginary values. So no divergences occur in the following 
models: 


• Timelike free bosonic theories: Momentum conservation dictates the value of 
the s-channel momentum. 


• S'L 2 (M) WZW model: In this model (reviewed in [1]), which describes strings 
in AdSs, large negative conformal dimensions require large values of the spectral 
flow number, which is not quite conserved but can take at most three values in the 
s-channel. 


So timelike theories can be hne, provided s-channel conformal dimensions are bounded 
from below in any given four-point function. This can occur if interactions impose re¬ 
strictions on s-channel quantum numbers. But there is no such restriction in timelike 
Liouville theory. 


3.3 Proposal for the spectrum and correlation fnnctions 

We now explain our proposal for the spectrum of Liouville theory with c < 1. In principle 
we keep the same real values for the momentum as for other values of c, but for computing 
four-point functions we move the s-channel momentum’s integration line slightly away 
from the real hne, in order to avoid all the poles. So the domain of integration in (3.5) is 

E = M. + ie , (3.9) 


for some real number e: 


4 - 

0 


Ps 

!i 


e 


The four-point function does not depend on the value of e. On the one hand, the four-point 
function does not change when the integration line is deformed in the upper half-plane, 
so long the behaviour at inhnity keeps convergent. On the other hand, negative 

values of e give the same result due to the invariance of the integrand under the reflection 

Ps —Ps- 

Then, instead of poles, the integrand of the four-point function has peaks: 
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Here we plotted the integrand as a fnnction of P = Pg — ie G W. The parameters are 
(3 = 1.103, (Pi) = (0.12,0.7,0.72,0.95),(;2 i,^2,^3,^4) = (0.4, 0 , cx), 1), and e = 3 x lO-^. 
Some of the positions of the poles are indicated in terms of the integers (m, n) snch that 
P = ^{m/3 + Some poles with small residnes do not prodnce noticeable peaks. 

Crossing symmetry can now be checked nnmerically. For example, with (Pj) = 
(0.32, 0.71, .45, .22) and {zi) = (0.27, 0, cx), 1), we hnd 


c 

s-channel 

t-channel 

precision 

0.873 

-3.68979 

-3.68975 

1 X 10-5 

0.564 

-3.17527 

-3.17526 

3.3 X 10-5 

0.241 

-2.74031 

-2.74032 

2.4 X 10-5 

-1.237 

-1.49777 

-1.49778 

1.1 X 10-5 

-3.751 

-0.50862 

-0.50861 

2.6 X 10-5 


Again, the precision can be improved by increasing the nnmerical accuracy of the com¬ 
putations. We can also study how the four-point function depends on the position x with 
(xj) = (x,0, 00 , 1 ). Taking (3 = 1.103 and (P*) = (0.18,-0.2,0.71,-0.43), we hnd the 
following plot: 
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We did relatively low precision calculations, in order for some difference between the s- 
and t-channel results to be visible. The difference is visible near x = 0 and x = 1, where 
the series representations of the t- and s-channel conformal blocks reach their respective 
radiuses of convergence. With our specihc choice of parameters, the difference happens 
to be more visible near x = 0 than near x = 1. 

Let us discuss the continuation of this four-point function to other values of c - that 
is, complex values of f3. As soon as (3 acquires a non-zero imaginary part, the poles of the 
integrand leave the real line, and inhnitely many poles cross our integration line: 


Ps 



- > 

0 



So although its integrand is analytic in (3, the four-point function (3.5) does not have an 
analytic continuation to complex values of f3. Actually, the existence of such an analytic 
continuation would have been an embarrassment of riches, as it would have provided 
a second dehnition of Liouville theory for c G C. On the other hand, the correlation 
functions are analytic in the momentums Pi, and can be computed for Pi & C and not 
just Pi G M. 

Finally, let us study the behaviour of the four-point function near zi = Z 2 - The 
universal behaviour of four-point conformal blocks is 

P^p\Pi\zi) ~ where 6 =-2P^ - 2P^. (3.11) 

^ Zi^Z2 2 
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On the other hand, P 2 , P 3 ) has a simple zero at Pi = 0 while C"^-^(Pi, P 2 , P 3 ) 

has a hnite limit. We must distinguish two cases: 

2 ;i 2 |'^| log| 2 ;i 2 |r 5 , (c<l) (3.12) 

2 :i 2 |'^| log 12 : 12!|~2 . (otherwise) (3.13) 

Therefore, the leading behaviour of the four-point function near zi = Z 2 depends on 
whether c < 1 or not. In both cases there are corrections in powers of Z 12 and 
The logarithms which appear here are due to the continuous spectrum, and have nothing 
to do with the logarithms of logarithmic conformal theory or of Coulomb gas correlation 
functions [14]. 

4 Non-analytic theories at rational values of c 

In our discussion of Liouville theory, we assumed that correlation functions are analytic 
functions of the central charge c and of the momentums of the helds. Under this as¬ 
sumption, and some other assumptions, there exists a unique non-rational conformal held 
theory at any complex value of c, which we proposed to call Liouville theory. 

We will now see that relaxing the analyticity assumption allows other theories to 
exist, when c takes the discrete values of eq. (1.1) - which we will call rational values. We 
give these theories the generic name of non-analytic Liouville theory, as their correlation 
functions are analytic neither in c nor in the momentums. The prototype of these theories 
is Runkel-Watts theory at c = 1 [15]. The generalization to other discrete values of c was 
proposed by McElgin [16]. 

We will check that McElgin’s proposal obeys crossing symmetry. We will also rule 
out other plausible proposals, including a puzzling family of three-point functions which 
obey crossing symmetry for some but not all values of the momentums. 

4.1 Non-analytic Liouville theory 

As we discussed in Section 2.3, the degenerate crossing symmetry equations have a unique 
solution provided b and b~^ are aligned and incommensurable complex numbers. If how¬ 
ever = - G Q, then there exist smooth functions which are periodic with both periods b 
and b~^. Using such functions, it is possible to build inhnitely many three-point structure 
constants which solve the degenerate crossing symmetry equations. Of course, we do not 
expect all these solutions to solve the full crossing symmetry equations. This nevertheless 
suggests that for ^ there may exist alternative solutions. 

An even more suggestive observation is that the DOZZ three-point structure con¬ 
stant, while it does not have a limit for generic values of c < 1 , does have a limit for 
rational values of c. This was hrst observed in the case c = 1 [3], where the limit of the 
DOZZ structure constant was found to agree with the already known structure constant 
of Runkel-Watts theory. This was then generalized to the other rational central charges 


X\Vp,{zi) 


li=l 


r\j 

Z\^Z2 


r\j 

Zl^Z2 


dPs 

dPs \Zl2\'^^^"P^ 
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[16]. The resulting three-point structure constant is 


^non-analytic^p^^ Ps, Ps) = Hm (Pi, P2, P3) = C"=^'(Pl, P2, P3MP1, P2, P3) , (4.1) 

^ p 

where we dehne, assuming that p and q are coprime integers, 


a(Pi,P2,P3) = 


_/l if n±±sin7r(i(p-g) +vW(Pi±P2±P3)) <0 , 


0 else . 


(4.2) 


This step function is the limit of , which is dehned for ^ M and invariant un¬ 
der shifts of momentums by b and b~^. We will now plot P 2 , P 3 ) and 

C"^-^(Pi, P 2 , P 3 ) as functions of Pi, and draw green vertical lines at Pi G where 

a{Pi, P2, P3) always vanishes. We take {p,q) = (7,5) and (P 2 ,P 3 ) = (0.2,0.45), and 
smooth the poles of into peaks by a small shift Pi ^ Pi + ie: 



We now consider the s-channel four-point function which corresponds to ( 7 n°°-anaiytic^ 


n 


i 2 = 1 


(P{—Ps,P3,P4) = l 


(4.3) 


We do not need to shift the integration contour away from the real line, as we did in Section 
3.3, because the poles (3.6) and (3.7) of the integrand all belong to |(/3Z-|-/3“^Z) = 2 ^=^) 
and are therefore outside the integration domain. This implies that our four-point function 
is a limit of the Liouville four-point function at complex central charges, as taking this 
limit does not involve poles crossing the integration domain. This in turn implies that 
non-analytic Liouville theory is crossing-symmetric. 

Near rational values of the central charge, Liouville correlation functions vary violently 
as functions of momentums. This makes it difficult to numerically check that their limits 
are indeed the correlation functions of non-analytic Liouville theory. So we focused on 
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directly checking that non-analytic Liouville theory is crossing-symmetric. For example, 
with {zi) = (0.23,0,cx),l) and (P^ = (0.22,0.37,0.28,0.12), we hnd 


c 

s-channel 

f-channel 

precision 

25 

2 

0.0245541 

0.0245559 

7.7 X 10-^ 

3 

5 

0.0547238 

0.0547235 

5.5 X 10-® 

0 

0.0361717 

0.0361713 

1.1 X 10-^ 

1 

2 

0.0909525 

0.0909523 

6.7 X 10-® 

1 

0.2937675 

0.2937677 

6.4 X 10-^ 


4.2 Other plausible theories 

Let ns relax the assnmption that p and q are coprime in the dehnition of the non-analytic 
factor a (4.2). For each rational valne of the central charge, we then obtain inhnitely many 
conjectnral three-point strnctnre constants, parametrized by pairs of integers (p, q) with 
a hxed ratio. Crossing symmetry is no longer gnaranteed, becanse the corresponding 
correlation fnnctions are no longer limits of Lionville correlation fnnctions at complex 
central charges. And we nnmerically hnd that crossing symmetry is not obeyed. 

Another idea is to replace the non-analytic factor a with any fnnction of cr, as was 
snggested in [17] in the case c = 1. Since a takes only two valnes, this amonnts to 
considering three-point fnnctions of the type 

a(Pi, P2, P3) = C"^-'(Pi, P2, P3) ((1 -y) + MPi, P2, P3)) , (4.5) 

for some parameter u, with Cq = and Ci = (^non-anaiytic^ s-channel decomposi¬ 
tion of the fonr-point fnnction reads 



' M+2e 


+ 1 /( 1 -z/) / + 1 / 


IEUE' 


' EnE' 


dPs 


C'^^'(Pl,P2,P,)C'‘=^l(-P„P3,P4) P^\p^\z,) , (4.6) 


which involves the disjoint nnion and intersection of the sets E and E' dehned by 


ls(P,) = a(Pi,P2,P,) , (4.7) 

1e'{Ps) = (t{-Ps.P:,,P^) . (4.8) 


As fnnctions of z/, the three terms of the fonr-point fnnction are independent polyno¬ 
mials, and we can investigate crossing symmetry term by term. The hrst and third 
term correspond to Lionville theory and non-analytic Lionville theory, and are therefore 
crossing-symmetric. It remains to investigate crossing symmetry of the second term - the 
integral over the disjoint nnion E U E'. Cnrionsly, we nnmerically hnd that the second 
term can be crossing symmetric or not, depending on the valnes of Pi, • • ■ P4. For exam¬ 
ple, taking (Pi,P 2 ,P 3 ) = (0.22,0.47,0.31) with c = 1 and {zi) = (0.23, 0, 00 ,1), we hnd 
the following behavionr of the s- and f-channel integrals as fnnctions of P4: 
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Based on this and other examples, we conjectnre that crossing symmetry is obeyed if 

((Pi.Ps) u (P3,pi)) n ((P2,pt) U (Pi.pi)) = 0 , with p, = ±s/pqPi mod § € (O, i) . 

( 4 , 9 ) 

This condition implies that the s-channel integration domain E U E' is made of intervals 
of the same lengths as the t-channel integration domain. 

In principle, crossing symmetry means the agreement of three decompositions, corre¬ 
sponding to the s, t and u channels. So far we have only considered the s and t channels, 
whose agreement for all values of the momentums implies that they also agree with the u 
channel. Now, in the case of Cip, we see that for any choice of the four momentums, two 
of the three channels agree. The third channel corresponds to using the operator product 
expansion of the two helds with the lowest values of Pi. This observation would deserve 
to be further checked, understood, and generalized to n-point functions. 

5 Conclusion 

Our results complete the dehnition of Liouville theory for any complex value of the central 
charge c. The spectrum is always given by the formula (2.9). On the other hand, the 
formula for the three-point structure constant depends on whether c < 1 (3.1) or not 
(2.27). 

Liouville theory can be characterized as the unique two-dimensional conformal held 
theory whose spectrum is continuous with multiplicities one, and whose correlation func¬ 
tions are smooth as functions of the central charge and momentums. Relaxing some of 
these conditions, we obtain other theories. In particular, we called non-analytic Liouville 
theory the theory which is obtained for rational values of c < 1 (1.1) as the limit of 
Liouville theory with complex c. 

We have also performed numerical checks of crossing symmetry in the generalized min¬ 
imal models, which exist at all values of c with the discrete spectrum S = 
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TZ(r,s)- In contrast to Liouville theory, generalized minimal models exist only on the sphere 
and not on higher genus Riemann surfaces, according to the following heuristic argument. 
The torus partition functions of generalized minimal models and Liouville theory are 
both inhnite. In Liouville theory, torus correlation functions become hnite as soon as 
non-degenerate primary helds are present. In generalized minimal models, all helds are 
degenerate and their torus correlation functions remain inhnite - except when they vanish 
due to the fusion rules. 

Now that we have a consistent non-rational conformal held theory for any value of 
c < 1, we can come back to our original motivation of accounting for statistical physics 
observables. Our results provide well-dehned predictions for correlation functions, which 
can be numerically computed to high accuracy for any values of the conformal dimensions 
and of the central charge. The challenge is to hnd matching observables in statistical 
systems. For example, it would be interesting to compare our results with four-spins 
correlation functions in the Q-states Potts model for general values of Q. Our results may 
also have applications in four-dimensional gauge theories, in the context of their relation 
with Liouville theory with a central charge less than one [18, 19]. 

A Numerical calculations 

In this Appendix we discuss which expressions for conformal blocks and correlation func¬ 
tions are suitable for numerical calculations, and which cutoffs and approximations are 
used. We do not specihcally describe our program for doing the calculations, but focus 
on the underlying principles and formulas. 

A.l Four-point functions 

An s-channel four-point function is an integral over an s-channel momentum Pg, see for 
example eq. (3.5) in the case of Liouville theory with c < 1. In Liouville theory, the 
integration domain is Pg E 'R + ie for c < 1, and Pg eR otherwise. In the latter case we 
can use the reflection symmetry Pg —)■ —Pg of the integrand, and halve the integration 
domain to Pg E (0, cxd ). 

We introduce a cutoff Pmax for the integral on Pg. Since the integral is exponentially 
convergent near Pg = oo, it is not necessary to take large values of Pmax- In the case 
c < 1 we have the additional parameter e, and the integral is over (—Pmax + fo, Pmax + fo)- 
Four-point functions should not depend on e at all. In practice e should be large enough 
for the integration line to be at a safe distance from the poles, and small enough for the 
integrals over the segments (Pmax, Pmax + fo) and (-Pmax, -Pmax + fo) to be negligible. 
We hnd that the four-point function is indeed constant over a wide range of values of e, 
typically 0(10-2) - 0(1). 

Our goal of computing correlation functions inhuences the way we compute structure 
constants and conformal blocks: 

• Conformal blocks are rational functions of Pg, with poles at the degenerate values 
Pg = P{r,s) for r,s E N*. We should compute the residues at these poles once and 
for all, in order to avoid recomputing them for each value of Pg when integrating 
over Pg. 
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• Structure constants depend neither on Zi, nor the conformal blocks’ cutoffs (such as 
the cutoff on (r, s) when summing over poles). We should avoid recomputing the 
structure constants when changing the values of Zi or these cutoffs. We however 
use built-in methods for integrating over Pg, and we do not control which values of 
Pg are used. This problem is solved by using spline interpolations of the structure 
constants as functions of Pg. 


A.2 Conformal blocks 

Global conformal symmetry determines the dependences of four-point blocks on three of 
the four positions Zi, and we assume (zi) = {x, 0, oo, 1). Let us write the conformal blocks 
in a way which is convenient for our numerical calculations, 


-r-(s')/r->l ^ _Sl_p2_p2/^ s_fil_p2_p2 

(Pj|x, 0, OO, I) = X 4 G ^2(^l—x) 4 G G 


N 


X 


( I + E E I . (A.l) 


7V=1 rs<N 


{r,s) 


where we introduced the elliptic variable q and function 6 * 3 ( 0 '), 

F{h i I, I - x) 


q = exp —vr- 




, Gsiq) = 


= 2^^ 

n£7j 


The coefficients are given by the recursion relation 


^(r,s) ~ -R(us) 


^N—rs,0 T ^ ^ 


^N—rs 


_ p2 _ p2 

r'"^N—rs , 


(A. 2 ) 


(A.3) 


where the prefactor is 


R 


{r,s) 


~2P(o,o)P(r,s) 


nr'=l-r ns'=l-s 2P(r',s') 


n n ± Pi + P{r',s')){P3 ± P 4 + -P(r',s')) • 

,2, it 
r —1—r s —1—s 


(A.4) 

In numerical calculations, we introduce a cutoff A^ax on the summation index N in eq. 
(A.l). We hnd that increasing A^ax by one often does not improve accuracy if Amax is 
even, so we take only even values of Amax- In this approximation, the non-trivial factor of 
the conformal block is a polynomial function of q of degree A^ax, and a rational function 
of P 2 with the poles {P^r,s)}rs<N^.^- 

Our formulas for conformal blocks are a reformulation of Ah Zamolodchikov’s recur¬ 
sion representation. (See [1] for a review.) Recursive formulas are particularly efficient 
for numerical computations, as compared to combinatorial formulas. The recursion rep¬ 
resentation which we used expresses a conformal block as series in q, whose radius of 
convergence is one. With the elliptic variable q, the duality which relates s-channel 
blocks with a position x and f-channel blocks with the position 1 — x has the self-dual 
value Q'seif-duai = 6 “’’’ ~ 0.0432. Since this value is small, there exists a regime where both 
s- and t-channel decompositions of the same four-point function involve conformal blocks 
with quickly convergent series expansions. 
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A.3 Structure constants 


The Liouville structure constant has different expressions (3-1) and (2.27), 

depending on whether c < 1. Both expressions are prodncts of T;, fnnctions, bnt the 
arguments of these fnnctions behave very differently. In the case of the argnments 
can have arbitrarily large real parts. So we can nse the integral formula (2.23) for the 
functions only after repeated use of the shift equations (2.22), and actnally it is simpler 
to use the product formula (2.24) instead. In the case of the argnments have real 

parts 'RQ or So we can more easily use the integral formula, which converges faster 

than the product formula. 

Let us start with the case of In the integrand of the fonr-point function, we 

have a factor Tfe(Q ± 2iPs), where Q ± 2iPs are at the edge of the domain where the 
integral formula for T;, converges. Using the shift eqnations, we find 

J] Tb(Q ± 2iP,) = 4P2 JJ r,{b ± 2iP,) . (A.5) 

± ± 


Using the integral formnla, we now find 


TT ^b{Q i 2iPs) 

y T6(f ±2*P,) 


4T6(6)2p2exp 



dt 4sin^ 2Psf sinh^ ^{b — ^)t 

t sinh bt sinh | 

0 


We will also use the identity 


(A.6) 


log 


nil T,(§ + 2tP,: 


n±.± Tj(2 + tPi ± iP2 ± iPs) 


= -4 


dt Ya=i sin'^ Pit - 2 sin^ Pit sin^ Pjt + 4 HLi 


t 


sinh bt sinh | 


, (A.7) 


which follows from trigonometric manipulations, using in particnlar 


sin^(ai ± 02 • • • ± an) = 2 


n—2 


cos 2ai 


(A.8) 


±,-± 


2=1 


The Ps-dependence of the prodnct of DOZZ structnre constants which appears in a fonr- 
point fnnction can therefore be written in terms of a single convergent integral, and we 
have 


^DOZZ(p^, p^)C^OZZ^_p^^ p^) ^ JJ 

2=1 


Tfe(Q + 2iPi) 
Tfe(f + 2iPi) 


xP^exp 

® Jot sinh bt sinh | ’ 


where we intro dnce the fnnction 


(A.9) 


ip{t) = sin^ 2Pst sinh^ |(6 — ■f)t — 2 sin"^ Pgt — (sin^ Pp — sin^ P2t)^ — (sin^ P^t — sin^ Pit)^ 
-I- sin^ Pgt i2 sin^ Pit — 4 sin^ Pp sin^ P2t — 4 sin^ P^t sin^ P4t | . (A. 10) 


2=1 
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If we instead used the integral formula for each T;, factor with a Ps-dependent argument, 
we would have to perform 10 integrals, which converge at t = 0 only thanks to regularizing 
terms in the integrands. 

In the case of the problem is the slow convergence of the infinite product (2.24), 
which is not very well approximated by hnite products with moderate numbers of factors. 
This can however be improved by approximating the infinitely many neglected factors in 
terms of an integral. To do this, we rewrite our inhnite product as 

CXD 

T6(f + iP) = exp ^ rij" , (A.ll) 

m,n=0 

where we introduce 

ym,n = ^ + mb + nb~^ , (A.12) 

and the function 

= log f 1 “ 4 • (^-^3) 

V y J y 

The idea is now to use the approximation 


CXD M N M pQQ N ^oo pco nco 

V~VV + V/ dn + v/ dm + / 

±^.JN+h t^nJM+h JM+l Jh 


dmdn . 


77T,,n=0 m=0 n=0 m=0 


n=0 -^^+2 


IM+IJn+1 


This can be made explicit using the following primitives of "0", 

f 1 \ 1 li % 

^'{y) = d log 1 + ^ + - + Pog-: , 

V y J y y-^ 


i’iy) = o + - 1) log 1 + ^ + iy\og 


2 2 


y^ 


y + i 
y - i 


which are such that '^(y) = 0 (-t)- We obtain the approximation 


y^oo 


(A.14) 


(A.16) 

(A.16) 


M N 


T,(f+ ,p) ~ A,-'” n If 1+ 


X e 


- 


m=0 n=0 


n=0 

ym.n 


y2 

am,n 




-Pb-iE)Lob' 


^M+^,n 


+ P=V' 


^M+^,JV+^ 


(A. 17) 


This converges acceptably fast, and allows us to accurately compute the values of using 
relatively small values of the cutoffs M and N. 
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